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jo. oralizn>d clescribin function is derived for a dissymetric 
relay with D» C, co. iponencs in both the input and output paths. Methods 
of evaluating tie D. C. components are discussed. The describing func- 
tion is applied to the analysis of step and disturbance inputs. 

The describing function is odified to allow computation of the 
D. C, and sinusoidal components of the error signal when a ramp input 
is applied. The applicability of the method in various situations is 
discussed. 

The validity of tne describing function analysis was investigated 
by comparing calculated responses with actual system response, as simu- 
lated on the analog computer. 

The authors wish to express their appreciation for the invaluable 
counsel and perceptive suggestions provided by Professor G, J. Thaler as 
faculty advisor for this project. 
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+ 


d 2 




2X 




X 
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d 1 


Pi 


+ 


P2 
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j:; i 

i::**TODucrio:. 

In receuu ./cars, tue aescribin., function ethod has been viidely 
used in the analysis of servo median is. is containing non-linear elements 
(1), The basic describing function technique consists of "linearising" 
the non-linear portion of the syste.i by assuming that: 

1 . The input to the non-linear device is a pure sine wave of 

known amplitude. 

2. The higher harmonics in the output of the non-linear device 

may be neglected. 

A representative describing function is then calculated by com- 
puting the amplitude gain and phase shift between the input and the 
fundamental frequency of the fourier series expansion of the output 
wave* These techniques are thoroughly discussed in standard texts (2). 

In some syste s, however, the presence of the non-linearity 
may cause a D. C. co lponent to exist in either the input or output 
path, or both. This is particularly true with dissymmetric non- 
linearities (3). A similar situation can be caused in many symmetric 
systems by introducing a ra ip input. In either case, the D. C. com- 
ponent may cause a change in the magnitude and/or phase of the funda- 
lental frequency of the output from the non-linear element, and it 
.may result in a significant change in the stability characteristics 
of the syste.i. 

It is therefore advantageous to define a generalised describing 
function, based on the following assumptions: 



1 



1. i.i_ . v. i : •!.- linear ele :eni. is a pure sine leave 

superi pose v " pure 2. C, component, 

2. i.ie higher har. onics in the output of the non-linear 

element :ay be ne glee ted. In this case, the output of 
the non-linear ele tent lay be reduced to a pure sine 
wave suneri. .nosed on a pure J. G. component, where 
botn ter :s .ay be evaluated from the fourier series 
expansion of the output wave. 

The generalized describing function still consists of the 
amplitude gain and phase shift between the sinusoidal component 
of the input and one f unda cental frequency/ component of the output 
wave. 

The purpose of this report is to develop the generalized de- 
scribing function for a dissymmetric relay and to investigate the 
characteristics of servo mechanisms containing such elements, fhe 
report includes a discussion of ;ethods of determining the magnitudes 
of the D, C. components; an analysis of system stability for step or 
disturbance inputs; and an analysis and discussion of the system 
response to step and ramp inputs. 
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* 'J.\ lx 

i. a.ig:: of d. c. sig al lsvsu 

2.1 General Considerations 

In all cases investigated in this report, it trill be assumed 
that the block diagram of the servomechanism can be reduced to the 
for : shown in figure 2.1, there G^s) and G 2 (s) represent the linear 
portions of a type 1 servo :echanism, and G^ is the linearized describing 
function which represents the non-linear element in the system. 




Figure 2.1 Block diagram of a servomechanism containing a dissymmetric 
non-linearity. 



In addition, it will be assumed that the servo loop acts suf- 
ficiently like a low pass filter so that only the first harmonic and 
the D. C. component of O(t) are transmitted around the loop. There- 
fore, for steady state operation with sinusoidal, step, ramp, or dis- 
turbance inputs, we can approximate the input to the non-linearity 
with equation 2. 1 ; and the effective portion of the output from the 
non-linearity can be approximated by either equation 2.2 or by equa- 
tions 2.3 through 2.5. 
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I(t) = '.-I + Xsin uit 



2.2 0(t) = -~ + a.cosujt + b.sinwt 

2 ‘ 1 

2.3 0(t) = ~ 4- Csin(tnt + <}>) 



2A 0 = J^ Z + b 1 2 



2,5 </> = tan - ^ ~ 

b 1 



It v:ill also be assumed that G^(c) represents a compensator 
having a transfer function of the form expressed in equation 2.6; and 
that G 2 (s) represents a rr.otor-load combination having a transfer func- 
tion of the form expressed in equation 2.7. 



Gi (s) 


= k 1 


7T(s + 


z i) 


8 




tF(TT 


Pi) 


where 


z i 


0 and 


Pi t 0 


G 2 (s) 


= k 2 


TT (s + z-i) 
1 s 77(s + p .) 


where 


z j 


0 and 


Pj / 0 



The basic Laplace transform equations for such a system ares 

2.3 £ = ©i - ©o 

2.9 i(t) =G 1 (s)(e i .e 0 ) 



2.10 0 0 = G 2 (s)0(t) 
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2 , 2 D, C. Le vel. x. i i-r Ut'-.'U' Path 

In the case of ste->, dis turbance , ran, or sinusoidal inputs, 
the D. C, co ioonent of the output fro : the non-linear element depends 
only on the linear pain of GpC 5 ) an d the applied input function, and 
is relatively easy to evaluate. 

The equation for the Laplace transform 0(t) can be obtained 
by using only the transfer functions of the linear portion of the 
system as in equation 2,11, Since the steady state sinusoidal com- 
ponents have no effect on the D, C, components in the linear portion of 
the system, the principle of superposition can be used to eliminate the 
transforms of these components when solving for the D, C. components 
only. The final value theorem can then be applied to obtain equation 
2.12, where 0(t)^, ©■ - , and I(t) c , represent the Laplace transforms of 
the non- sinusoidal components of 0(t), , and I(t), respectively. 



2.11 



0(t) 



_ _ 3 0 - 

u^IsT 



VCsT 



6 i 




2,12 




s=o 



^2< s ) 




I(t) d 
G-| ( s ) 



s=o 



In order to simplify the following calculations, the limits of 
G^(s) and G^s) as s approaches zero can be partially evaluated to 
get equations 2.13 and 2.1':. These results can then be subs Mb it ad 
into equation 2,12 to get equation 2,15. 



2.13 G^s) 



S-*-0 



= v 1 TT(s * Z i } 
77 (s p i ) 



s-*o 



_ , 77 z 1 _ ^ 

‘TTVi S 



5 



/ \ 

\ ^ / 



Z *> O 



77 (c 
077 



S*-C 



77 s.- 



77 



-o - 
2 " 



s 2 
r — 
“ v 2 



© 



id 



I(t), 



V 1 



S.f o 



for a step input, or disturbance, of magnitude a, the Laplace 
transfer * of is given oy equation 2,1, For a steady state analysis, 
the initial tine is in ate rial, and consequently the Laplace transform 
of I(t) d can be represented by equation 2.1? and the D. C. component of 
C(t) can therefore be computed by equation 2.1". 



2.16 0 . = Au(t) 


_ A 




1 


s 




2 . 1 ? Kt) d = j 






2 in a o _ c 2 

a, le ~ - rr*— 


" " 


= 0 


2 




s^-o 



a similar argument can be applied to sinusoidal inputs, since 
in this case, &id is zero. 

For ramp inputs of magnitude 3, the Laplace transform of the 
input is given by equation 2.19, and the D. C. component of 0(t) can 
be evaluated by equation 2,20. 

2.19 e i = 3tu(t) = - 2 - 

s d 
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For an n° n order 



he s ler. in the denominator of 



Gp ( s ) would be changed to s 1 Jquacion 2,21 can be derived for such 
a system by following a line of reasoning similar to that used for the 
first order system, Evaluation of equation 2,21 shows that the d.c, 
component of 0(t) will be zero for both step and ramp inputs for any 
system of order higher than one. 



2.3 D« C. Levels in the Inout Path 

In the preceding discussion, there was no limitation on the 
nature of the non-linearity represented by G^, since the D, G, com- 
ponent in the output path is independent of G^ , In general, however, 
the D. C, level in the input path will not only depend on the nature 
of the system input and the characteristics of the linear portion of 
the system, but will also depend on the nature of the non-linear ele- 
ment , 

In many cases the solution for the D, C, component in the input 
path can be extremely complicated, and nay result in transcendental 
equations which can best be solved by graphical means (3), 

It should be noted, however, that in many cases it is unnecessary 
to solve explicitly for the D, C, level in the input path in order to 
obtain the describing function curve. This fact is particularly well 
illustrated in the analysis of the relay servo with a ramp input, given 



2,21 



a o _ s n ft 

T ~ It 




in section 4,2, In other cases, such as the analyses for step inputs 
given sections 3.2 and 3.3i the solution for the D. C. level indicates 
that for frequency responses analysis, the D, C. level in the input path 



ca.. us ±3* in' li’i !»■'*. • r-i »/ is treated as a cy 

:etrical co pontirc, 



c:Ar j:i in 

TIU RELAY 3ER'/C 

3.1 Discussion of the Describin'- Function . 

Appendix I contains the derivation of the generalized describing 
function for a dissy 'metric relay with the following characteristics: 

1 . l’he pull-in voltage for a positive input is not equal to 

the pull-in voltage for a negative input. 

2. The drop-out voltage for a positive input is not equal to 

the drop-out voltage for a negative input. 

3. The hysteresis for a positive input is not necessarily equal 

to the hysteresis for a negative input. 

4. The positive relay output is not equal to the negative re- 

lay output. 

The input vs. output characteristics of such a relay are shown 
in figure 1-1 of Appendix I. 

Since most practical relays have positive and negative outputs 
of equal magnitudes, no attempt was made to analyze a relay with un- 
equal outputs. Therefore, in the remainder of this report, it will be 
assumed that the relay output voltages are of equal magnitude. 

3.2 The Describing; Function for a Relay vrith no Hysteresis . 

The input vs. output characteristics for a dissymmetric relay 
with no hysteresis can be represented by figure 3.1. 

Equation 2.13 shows that for a servomechanism of the type shown 
in figure 2.1, the D. C. component of the relay output is zero for 
steady state operation after the application of a sinusoidal, step, or 
disturbance input. Combining this result with equation 1-15 of Appendix 
I gives: 



p 




Figure 3. 1 Transfer characteristics of a dissymmetric relay with no 
hysteresis. 
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.i 

TT ' 



^ 1 - &2 



The angle relationships of equations 1-7 and 1-11 of Appendix I 
can hot: be introduced to solve for the D, C. co lponent of the input 
wave and to solve for A and /?•-> as shown below. 

0. . — W n do + W 

3.3 smftj = _!_ = sxnR = — 



3.^ 



d i_ d 2 

2 



3.5 sin^j = = s p _ # 



r; 



2 



For a relay with no hysteresis the expression for a., as derived 



in Appendix I, is 



3.5 a.j = 0 



hie expression for b- is given by equation I-lo of Appendix I. 
co Ibining this e qua tier with equation 3.2 gives 



3.7 t. = ^.(cos B 9 -!- COS $ ) = z±cos£ 

7T 1 77 1 



4Y 

nr 



d 1 J d 2 



2X 



/' - S -* B \ =? i 1 - 

the equation for the generalized describing function is; 

-1 



3.3 G, = I / a " ; b 



2 + b 2 

1 °1 



lan p. 





For the dissymmetric relay with no hysteresis, the describing 
function is; 
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for 



> 



Oh 



d 2 



9 



Table I is a list of the values of the generalized describing 
function for a relay vjix.li d^ + d£ = 5.2, and ! = 1 . A graph of this 
describing function is shown in figure 3.2. 



X 


n 


G d _1 


Of 1 (db.) 


2.6 


0 


OO 


o° 


3.0 


.212 


4.73 


13.5 


3.6? 


.257 


3.89 


11.5 


4.0 


.242 


4.13 


12.25 


5.0 


.217 


4.60 


13.2 


6,0 


.191 


5.23 


14.3 


7.0 


.169 


5.92 


15.4 


8.0 


.150 


j 6.65 


16,4 


10.0 | 


.123 


S.13 

1 


17.2 



Table I Values of Gj for a relay with no hysteresis, 
d 1 + d 2 = 6.2 

Equation 3.9 applies only for cases in which the relay operates 
on both sides during each cycle. It is not necessary to consider one- 
sided operation at this point, however, since this type of operation 
requires that the 3. C. level in the relay output .lust be different from 
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zero. Io has -reviousl,, ^ee.. shown Chat for a servo .icchanis 1 of the 
type described in Chapter II, the stead state D, C, component in the 
relay output vast be aero for sinusoidal, step, or disturbance inputs. 
If the D. C, level in the output is allowed to be different from 
zero, as in the case of a ramp input, then both one-sided and two- 
sided operation is possible, The subject of ramp inputs is treated in 
Chapter IV, 

3.3 The Describing Function for a Relay with hysteresis . 

The input vs, output characteristics for a dissymmetric relay 
with hysteresis can be represented by figure 3.3 • In this case, the 
combined results of equation 2,13 from Section 2,2, and equation 1-12 
from Appendix I gives: 



3.10 



a o - 2L(* 
”2 2 IT 2 





3.11 /£+<,=^ 2 + 



Combining the angle relationships given by equations 1-3 through 
1-11 with equation 1-13 of Appendix I, we obtain equation 3.12, 



3.12 &1 = ~ 



d 1 " Pi + y2 " P2 



A 



If the hysteresis is the same for both sides of the relay, 
d.j - p.j = d 0 - p, , and the angle relationships can be evaluated as 
follows : 

3.13 sin A - sin c*'. = Hi = H.2 ’2 

1 1 X X 

= sinfi - sino^, 
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Figure 3 . 




lot * «(, 



to t 3 IT - 






lot «2fT-^j 



3 Transfer characteristics of a dissymmetric relay with 
hysteresis. 



. 1 sin^j - si: ^ 



2sin:;(^ - o^cos-K^ + <j) 

2sin-K^ 2 " ^2^ COS ^^2 + ^ 2 ) 



By combining the results of equations 3 • 1 1 » 3-13* and 3- 1 if ", it 
can be shown that ^ ^ and 8 ^ - /? 9 < The D. C, component of the 

input can then be evaluated by equation 3.l6» and the angles can be 
evaluated by equations 3*17 and 3. IS. 



5.15 sin/ 9 . = — _ _ 

A 



a 2 



= sin 4 



3.16 ■ = d 1 - d 2 ?1 - ?2 

2 2 



3.17 sin^ = sino< 2 = 



. 1 sin^j = siu^, = ~1 



- Pi * P2 
2X 

di -!- d 



= R„ 



2 _ 



2X 



R, 



The general expression for b^ is given by equation 1-1 4 of 
Appendix: I. Evaluating this equation by inserting the results from 
equations 3.17 and 3. IS gives equation 3.19. The describing function 
equation, 3.20, is then obtained by substituting equations 3.12 and 
3.19 into equation 3.8. The describing function in this instance also 
reduces to exactly the same form as the describing function for the cor- 
responding symmetric relay. 



3.19 b. = — (cos^i + cosO 

1 77 1 1 



= f?< J 1 - % 2 + 7 1 - -v 2 > 



16 




4 > 









i - r 3 2 



for X > ?1 * ?2 
2 



Table II contains the values of the generalised describing func- 
tion for a relay with + dg = 6 and pg = 10. A graph of this 

describing function is shown in figure 3.^. 



1 ? 



■< r 

A 


G d 


V 1 


\i~ ^ (db) 


Phase Angle 


5 


.1135 


» 


18.9 


26,6 


5.2 


.1413 


7.03 


17.0 


19.5 


5.5 


.1505 


6.65 


16.4 


16.2 


6 


. 1 54: 


6. 46 


16.2 


13.2 


7 


.143 


5.75 


16.7 


10.1 


n 


.1372 


7.28 


17.2 


8.33 


9 


.1265 


7.9 


17,9 


7.15 


10 


.1163 


r'. 

^00 


13.7 


6.26 


11 


.103 


9.27 


19.3 


5.62 


12 


.100 


10.0 


20 


5.1 


15 


.03018 


12.22 


21.7 


3.97 


20 


. 0624 


16.0 


24.1 


2.92 


30 


.0422 


23,7 


27.4 


1.92 


40 


.0316 


31.6 


29.9 


1.44 


50 


,025-4 


39.35 


31.8 


1.15 


Table II 


Values of 


G^ for a relay with 


hysteresis. 






d 1 + d 2 = 


6 and P-] + P 2 = 10 
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(S& n& f ajiz. eJ 



Function 

** G K 



£& ser/jo inf 

Uj! 

,-t-t 



25 



v.io sa.ie for booh sides of the relay, 
ohe uranscei. eiital e .uatishs ohici .ust be solved in order to determine 
the D. C, co ponent of I(t) become m,ch more complex. They are difficult 
to solve with Lrigono 'etric relationships alone, but it is relatively 
easy to solve then by graphical . ehhods using the construction shown in 
figure 3,6. This construction can be developed in the following manner: 

1 , Selected a prooractor with radius ?.. 

2, Calculate the following quantities for each value of X to be 

investigated. 



= ’-’I' 



- <-1 



>2 



- d 2 ?t 



O 

J* 



Construct figure 3.5 using the calculated dimensions. 

4, hove the protractor along the vertical axis until 

$1 + ^ = A + as shown on figure 3.5. This satisfies 
equation 3.11. 

The construction guarantees that equations 1-3 through 1-11 are 
satisfied if and only ^ and are selected as shown in figure 
3.6, and the D, C, component of I(t) is established by the center of the 
protractor, figure 3.6 therefore represents the desired solution for the 
angles, and b = 

The generalized describing function can then be evaluated using 
equations 3.3, 3.12, and 1-14. Table III gives the values for such a 
describing function for a relay with p^ - 3, Vo ~ 7, d^ = 2 and = 4, 
xhis desci’ibing function is plotted in figure 3.7. 
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Figure, 3.5 




Figures 3.5 and 3.6 Graphical construction for solution of equations 

3.11 and 1-8' through' 1-1 1 . 
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4 


G d “ 1 (db) 


.2 


-1. 00 


, 12°0 


21 .4 


17.3 


6.0 


-1.615 


.1575 


13.45 


16.25 


7.0 


-1.p62 


.1525 


10.2 


16.25 


8.0 


-1 . y-0 


.1370 


co 


17.1 


10.0 


-1.520 


,1164 


6.3 


13.? 


12.0 


-1.520 


.0999 


5.1 


20.0 


15.0 


-1.510 


.0813 


4.0 


21.7 


bote : 


Values for and 


P for larger 


values of X are 


approximately 



the same as corresponding values in table II, 



Table III Values of the generalized relay describing function for 
P-j = 3* P 2 = 7, = 2, and d 2 = 4. 

3.4 Cxmarlson of Sy metric and Dissymmetric Relays . 

Equations 3.3 and 3.20 indicate that the generalized describing 
functions for dissy metric relays are exactly the sane as the corresponding 
describing functions for syr .etric relays, provided the dissymmetric relays 
have equal hysteresis on each side. A s e rvone chanisn containing such a 
dissymmetric relay would therefore respond to a sinusoidal, step, or 
disturbance input in a manner quite similar to a servo containing a sym- 
metric relay. 

If the servo response results in a limit cycle or sinusoidal out- 
put, the output oscillations for the symmetrical relay would be centered 
about the position commanded by the D. C. portion of while the output 
oscillations for the dissy metrical relay would be offset slightly to 
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[GENCRALIZ 

DESCR I3ING F 



UNCTION 



DISSYMMETRIC 



X »-/2 

l i : 

\Xs/0 



create the .ncessar u, Z t me". ach signal at the input of the relay, 
rhus the effect of hiss;* e try is to introduce a steady state D, C. error. 
If the diss/ metric relay does not have equal hysteresis zones on 
each side, the generalized describing function nay be quite different 
fro:: the describing function which would be obtained by assuming equal 
hysteresis zones and using equation 3.20. A comparison of figures 3.4 
and 3.7 indicates that these differences will be relatively minor in most 
practical cases, and that adequate results can be obtained by using equa- 
tion 3,20. 

3.5 Verification of Proposed Describing F-.mctior.s . 

In order to verify the describing functions previously derived, a 
position control system consisting of the dissymmetric relay and third 
order linear components was postulated. Linear transfer functions with 
typically low pass characteristics were selected for this system. 

Figure 3*3 is a schematic of the control system together with the 
analog computer simulation. This simulation provides for operation of 
the relay under varying degrees of dissymmetry and hysteresis. (4) 

Figure 3*9 is a bichols Plot of the linear transfer functions with an 
arbitrary gain factor of 20, 

Oase I Dissymmetric Relay Without Hysteresis. 

The relay xias operated without hysteresis and with a constant 
total dead zone of 5.1 volts. The plus and minus pull-in voltages were 
varied on a series of runs to obtain various conditions of dissymmetry, 
3tep inputs of 10 volts and 15 volts were applied to the system. The 
gain for the linear portion of the system was adjusted to different 
values so as to place the linear curve in the desired position on the 
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TmJfrWi 




.iciiols hlo 



w, Ui ,-o .v i. a series of lfLiit cycles of various 

amplitudes , 

fable IV is a siuriary of tie results obtained, Typical brush re- 
cordings fox- several values of system gain and input magnitudes are dis- 
played in figures 3.10 through 3.12, which illustrate the predicted steady 
state d.c, error levels and sinusoidal characteristics. 



Run 


Relay 

Pull-in 


3yste:.i 

Gain 


Observed 

Amplitude 


Li lit Cycle 
Frequency 


Predicted 

Amplitude 


Limit Cycle 
Frequency 


1 


+2. 6, -2. 5 


QO 

• >o 


4.6 


2 


4.7 


2 


2 


+3. 1 ,-2,0 




4,6 


2 


4.7 


2 


3 


+2.4.-2.7 


132 


6.0 


O 

/A. 


3.0 


2 


4 


+4.0.-1.1 


132 


7.9 


2 


0.0 


2 


5 


+2.0.-3.1 


132 


3.1 


2 


0,0 


2 


6 


+3.5.-1.6 


132 


8.0 


2 


o 

CO 


2 


7 


+2. 7, -2,4 


81 


4,0 


2 


3.3 


n 


8 


+2, 7, -2.4 


63 


-.one 


- 


lone 


- 


Table 


17. Experimental 


Results for 


Dissymmetric 


Relay Tithout hysteresis. 



Case II Dissymmetric Relay with Rystei’esis. 

The relay was overate* 1 with the following characteristics: pull 

in, + .93 volts and - 3. volts; drop out, + 4.0 volts and - 1.93 volts. 
The sane general procedure was followed as in Case I. A sunnaiy of results 
is shown in table V. typical rush recordings are displayed in figures 
3.13 through 3.13. 

A co iparison between the predicted and observed values shown in 
table 7 indicates the validity of the proposed describing function for a 
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relay wi cl hysteresis. It n .... . e ,,o oeC. that a small difference between 
the observed and predicted values is L-o be expected since the experimental 
relay characteristics wore nod exactly matched to those used for predic- 
tions. file predicted values mere based on figures 3.^ where p-j + pg = 10 
and d-j + dg = 6. the characteristics of the experimental relay were : 

Pi + p 2 - 9.2 and d^ + d 9 = 5.95. 



Run 


Gain 


Observed 
A rplitude 


Limit Cycle 
frequency 


Predicted 

Amplitude 


Limit Gycl* 
Frequency 


1 


620 


39.6 


1 . '3 


*}0 


1.93 


2 




29.0 


1.38 


30 


1.90 


3 


300 


19.35 


1.75 


20 


1.35 


k 


1?2 


12,75 


1.72 


12 


1,75 


5 


3? 


7.15 


1.50 


6 


1.51 


fable 


V. d: 


■•perl lental Results 


for Dissyrr 


metric Relay 'Jith 


Hysteresis . 
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5 of- st-erj responses. 
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4. 1 General Gonsidera .dons . 

..lien a relay servo is subjected to a step or disturbance input, 
the servo '..'ill either eo a oo rest while the error remains in. thin the 
lriits of the dead zone, or it will enter a limit cycle which operates 
each side of the relay during each cycle, binilarly, when a sine input 
is applied, the steady state error signal will either vary id thin the 
limits of the dead zone, or it will cause the relay to operate on both 
sides during each cycle, 

’.ihen subjected co a ra ip input, however, the error wall never re- 
rain in the dead zone, the servo output .ay core to rest, but the in- 
creasing ramp signal in uhe input channel '..dll always cause the relay 
to actuate again. the ref ore, the servo rust enter a Unit cycle when- 
ever a ra:rp input is applied, unless the ramp signal is so largo that 
the relay does not drop out. this limit cycle nay involve actuation 
of both sides of the relay, or it ay involve only the repeated operation 
of one side. 

The case in which the relay operates on only one side is the 
simplest, and will be discussed in lection 4.2. hen the relay operates 
on both sides, the solution must be obtained by graphical methods, as 
explained in section 4,3. 

4. 2 Describin'; Function for One-sided Operation . 

.hen the relay operates on only one side, the input vs. output 
characteristics can be jorurayed by figure 4.1, The fourier expansion of 
the output wave, ‘which is derived in Appendix I, becomes 
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Figure 4.1 




Transfer characteristics of relay operating on only one side. 



\ 



3 ? 



.1 






7 T 



TT - C°c /3 ) 



' a-1 - £ (si i£- sino() 

! 7T 



' .J = " (cos /9 + COS oO 



.Siei. the servo has a ramp input, ©. = btu(t), the D. C. component 
of the relap output irave is pi von bp equation 2,20. fhis equation can 
be co ibined T .rith equation 4. 1 to pet equation 4.4 



2.20 fft = 3 

^ v 2 

4.4 = Aor-*-*) 

-^V2 ^ 

The equations for the generalized describing function can now be 
obtained in the following manner, without solving explicitly for the D. 
C. component of the relay input. 



4.5 a^ 2 + b^ 2 - 



TT 

pi 2 

TT 



= 2 



= 2 



2 

’'(sin 2 /? + sin 2 o< - 2sin$sina0 
(cos £ + cos o( + 2cosificoso() 

(2 + 2coso(cosB - 2sin o( sin B ) 

2 



1 + cos(<\ + 



1 - cos 



1 - cos 



(7 T- <£ -yfl)J 

22 n 1 

•vy I 

“Vo J 



4.6 tl = 52bi - sliLgC - ^ 1( / . ^ ) 

b^ cos/9 + cos 2 

4.7 ^ = tan - 1 1± A(/S-*) 
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ais . , ii xs convenient to introduce the trigonometric 
identity given by equation 4.8. 



4,3 sin£~ sinc( - 2sin i(^-o()cos + o() 

2 2 



= 2 sin p cos 



<7r 

2 



bit 



in. 



■V2 



= 2 sin <j> sin 



YK- 



■v- 



4.9 sin p = SniL=_£Hl* 



2sin -. w 
Xhv 2 



4.10 sin / = 



d i - ?i 
2Xsta J2. 
" v 2 



The generalized describing function for the ramp input during 
one-sided ooeration becomes 



' s - 11 0 d‘5f* 



1 - cos 



23 7T 



YIC 



v 2 



0 



The relay will operate on only one side as long as X is less 
than a critical value, X c , defined by equation 4.12. 

b 12 r - Pi + p 2 

4.12 X c 



1 + sin oi. 



If the relay has no hysteresis, there should be no discontinuity 
in the describing function as the servo shifts from one-sided operation 
to two-sided operation. On the other hand, if the relay has a signifi- 
cant amount of hysteresis, there will, in general, be a discontinuity 
when the shift occurs. 
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contains no .niniaum li.nl o 



The ex. r«ssion lor u .e magnitude of 3^ 
for the amplitude of the sinusoidal input* The phase angle expression, 
however, indicates that the describing function is valid only if 



L'ote that if the denoininator of this expression is equal to its maximum 
value, 2, this restriction lerely requires that the peak to peak oscil- 
lations of the sine wave must extend from one side of the hysteresis zone 
to the other. 

For relays with absolutely no hysteresis, the describing function 
predicts absolutely no phase shift due to the relay. However, if the 
relay contains at least an infinitesimal amount of hysteresis, as all 
real relays do, there will be some small value of X for which the magni- 
tude of G^ is defined, and for which the phase shift will be 90 degrees. 
Consequently, the effect of hysteresis must be considered to be much 
more significant in the case of ramp inputs than for other types of inputs. 

The describing function defined by equation 4.11 can be plotted 
in non-dimens ionalized form by defining the following quantities 





4.14 





4.15 sin f 



23sin _JT 
—-v 2 



I— Cm(s) m(s) 

1 " d 1 



' *1* O o 

4.1? D ? , - _J 1 



The critical equation for stability becomes 
4.13 3- d :. f1 (s)l 2 (s) = 4 dn l(s) = -1 



The plot of 



J d> ; ! db 



hZl 



can then be used as che ra.np describing 



function for one-sided operation of any relay by simply adjusting the 



gain of the linear plot by the factor 



? 1 - a. 



. A graph of the de- 



scribing function of equation 4.14 is plotted in figure 4.2. 

.'he naximun value of 2 for which the describing function will be 
valid is 

4 10 7. = , :: c = Pi + ?2 

, y c ppnq- T^”7q7Ti“+-sTH^a' 

— - ^hc 

1 + sin 

The describing function plot can also contain loci of equal 
values of = Z(1 + sin o { ) . These loci will pass through points which 
correspond to the naxi an value of 2 for which one-sided operation trill 
occur for a relay with the corresponding ratio of dead zone to hysteresis. 

In order to ootain these loci, it is rost convenient to first plot 



Z vs Z(1 + sin o() for various ratios of 



This requires solvin: 



‘ v 2 



for the angle $ using equation 4.15 then solving for the angle o( using 
equations 4.4 and 4.7. The values of Z can then be obtained for constant 
values of and these points cam then be plotted on the describing function 
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i 

t= 



Figure 4.5 Illustration of wave shape of I(t) in- non-sinusoidal region 
of ranp describing function plot. 
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cur .e 



iS 



values of 3 . 

1 Ai/ 2 

curve for 3, = 10, 



U'c , IolS o_ 3 vs «(1 + sin o ( ) for various 
, inf or .atio’i on these curves was used to obtain the 
xjhich is sno’.jn on figure ’+.2. 



It can also be shown that there is a li- tilting value of phase 



shift, beyond which the relay input can no longer be sinusoidal. If 

3 > , 5YIC , sinusoidal operation is possible only if 
“2 



4.20 -j£ S 77- J£L 

^2 

If 3 ^ .51 K , sinusoidal operation is possible only if 
v 2 



4.21 -£• < «j!L 

v 2 

These Uniting values are also plotted on figure 4.2. The 'mathematical 

proof of the limitation is contained in Appendix II. 

The non-sinusoidal characteristics become more and more pronounced 

as — approaches zero or unity, and as $ approaches 90 degrees. For 
m V2 

example, for 3 = .1 Yk v and j£ = 81°, the solution of equations 4.4 and 

2 

4.7 requires that ck - 1?2 C and & - Q° , This requires a wave shape of 
the form illustrated in figure 4.5. The corresponding sinusoidal wave is 
shown in dotted lines for comparison. 

4.3 Experimental Procedures 

As a first step toward the determination of a describing function 

for ramp inputs, a series of investigative runs were made on analog 

control systems. For this purpose two control systems were simulated with 

transfer functions ^ . , and , " , respectively, hichols 

3X8 + 4) .J(X + 2)(.b + 4; 



Plots for these systems are shown in figure 4.6. 
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*inaipsis o. w/.e results of theso preliminary runs together with 

the theoretical concepts uhich had been developed up to this point, 

indicated the significance of a oaraneter formed by the ratic — , 

111 ,, 



An experiment was then designed to ’lore completely explore this feature, 

A series of values for the ratio — — , ranging from 0.01 to 0.90, 

Xa V2 



were selected. Several runs were made at each of the selected values 
of the ratio, by varying B and K to keep the ratio constant. The 
relay was operated at a constant value of hysteresis and data was 
recorded for both the second and third order systems. The reduced data 
is shown in table VI, and on a Nichols Plot in figure 4.7. 

An examination of the data plotted in figure 4.7 leads to the 
formulation of certain conclusions. Since each set of limit cycle condi- 
tions, plotted for a constant value of % s seems to lie on a single 

rKv 2 

curve it appears that a describing function curve exists for each value 
of this ratio. The fact that the limit cycle conditions obtained for 
the second order system substantiate the pattern established by the 
third order system even more clearly illustrates the existence of a 
describing function. The nature of these describing function curves in- 
dicates that they will extend across all phase angles from -90° to -180° 
if the curves are extrapolated beyond the region of plotted data. This 
agrees with intuitive reasoning that a limit cycle must exist for all 
type I systems with ramp inputs, unless the system is incapable of fol- 
lowing the ramp. 
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X 
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rad/ sec 
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db 
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volts 
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db 


.5 


.01 










.924 


2.0 


103 


34.5 


.5 


.02 


.342 


4.45 


109 


36.2 


1.08 


1.47 


106 


26.7 


1.0 


.02 










1.515 


1.94 


111 


29.9 


.5 


.05 


.533 


2.40 


125 


23.5 


1.33 


0.95 


109 


16.5 


1.0 


.05 


.61 


3.67 


130 


28.2 


1.965 


1.21 


116 


19.0 


1.5 


.05 


.633 


4.8 


134- 


30.5 


2.40 


1.4 


121 


20.2 


2.0 


.05 










2.73 


1.56 


124 


21.3 


2.5 


.05 










3.14 


1.67 


128 


21.8 


.5 


.1 


.64 


1.65 


131 


15.7 


1.65 


.72 


112 


8.9 


1.0 


.1 


.796 


2.26 


140 


19.1 


2.325 


.86 


120 


10.9 


1.5 


.1 


.862 


2.95 


142 


21.7 


2.86 


.97 


125 


12.2 


2.0 


.1 


.932 


3.51 


146 


23.3 


d.p 1 


1.10 


129 


13.1 


2.5 


.1 


.96? 


4.70 


147 


24.6 


3.59 


1.29 


132 


14.1 


3.0 


.1 


1.03 


4.60 


149 


30.1 


3.92 


1.22 


134 


19.1 


1.5 


.3 


1.16 


1.47 


155 


8.3 


3.31 


.62 


129 


1.2 


3.0 


.3 


1.45 


2.10 


166 


10.5 


4.72 


.74 


140 


2.3 


1.5 


.5 


1.12 


1.08 


153 


4.0 


2.86 


.57 


125 


-1.3 


2.0 


.5 


1.21 


1.26 


156 


5.6 


3.40 


.55 


130 


-1.3 


2.5 


.5 


1.35 


1.43 


161 


5.8 


3.99 


.60 


135 


-1.1 


2.0 


.7 


.932 


1.11 


147 


5.7 


2.415 


.57 


121 


- .1 


2.0 


.9 


.472 


1.16 


122 


11.7 


.967 


,62 


103 


6.6 



Table VI. Lata obtained in experimental determination of describing 
function for ranp input. See figure 4.7 for plot. 
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Another 



L-restin^ fact observed from the plots in figure 4.7 is 



cvrnrastry ax>ut the value = 0,5. Thus the curve for = 0.4 

o o -tr-r- v t ' 

XaV- Y i\-_ 



YI 



' v 2 



coincides with that for 0,6, the curve for 0,3 coincides with that for 
0.7 etc, 

4.4 Verification of the Proposed Describing Function 

Although the experimental results presented in section 4.3 con- 
firmed the existence and to some degree the nature of the describing 
function, another experiment was designed to explicitly verify the 
theoretically derived describing function. In addition to the second 
and third order control systems previously simulated, a, fourth order 



system with transfer function 



K 



was used. In con- 



S(S + 1)(s + 2) (S + 4) 
ducting these runs the gain of each linear system was held constant so 

that the linear curves could be superimposed on the plot of describing 

function curves, and the predicted limit cycles were obtained from the 

intersections. The value of the ramp input was varied in order to obtain 



different values of the ratio JX 



YIC 



The hysteresis was set at approxi- 



’ v 2 



mately one volt and held constant throughout. Yichols Plots of the 
three linear systems are shown in figure 4,6 and the nondimensionalized 
describing function curves are shown in figure 4.2. The predicted and 
observed limit cycles for each system are compared in table 711. The 
Brush recordings are displayed in figures 4.9 through 4.21. 

Since some of the predicted limit cycle values showed dis- 
crepancies of 25 . j or more from those observed, it is felt that a dis- 
cussion of error sources is warranted here. First of all the inherent 
errors of any describing function analysis are present, that is the 
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X 
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Freq, 


1 


.05 


.281 


3.16 


.209 


3.29 


.150 


3.9 


23.9 


2 


.10 


.562 


3.00 


.366 


2.13 


.44 


27.3 


20.2 


3 


.20 


1.125 


2.89 


.61 


2.12 


.75 


25.8 


23.0 


4 


.30 


1.687 


2.79 


.748 


2.35 


.07 


15.7 


8.2 


5 


.50 


2.812 


2.73 


.350 


2.56 


.93 


7.6 


9.4 


6 


.70 


3.937 


2.75 


.722 


2.35 


.37 


14.5 


21.0 


7 


.80 


4.400 


2.85 


.572 


2.12 


.75 


22.0 


30.0 
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Order 
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1 


.05 


.231 


1.74 


.379 


3.21 


.170 


84.0 


55.0 


2 


.10 


.562 


1.60 


.679 


1.77 


.56 


10.6 


17.5 


3 


.20 


1.125 


1 .49 


1.17 


1.32 


1.22 


11.4 


4.3 


4 


.30 


1.687 


1.4-0 


1.48 


1.32 


1.60 


8.6 


8.1 


5 


.50 


2.812 


1.36 


I .67 


1 .40 


1.35 


2.9 


10.4 


6 


.70 


3.937 


1.45 


1.39 


1.32 


1.60 


9.0 


14.3 


7 


.80 


4.400 


1.53 


1.07 


1.32 


1.22 


13.7 


14.0 








Second Order System 








1 


.05 


1.0 


1.21 


1.965 


2.52 


.75 


108 


62.0 


2 


.10 


2.0 


1.10 


3.31 


1.56 


2.4 


41.8 


27.5 


3 


.30 


1.5 


.62 


3.31 


.580 


3.37 


6,4 


1.8 


4 


.50 


1.5 


.57 


2.86 


.490 


3.00 


14.0 


4.9 


5 


.70 


2.0 


.57 


2.41 


.504 


2.40 


11.5 


0.4 


< 

O 


.90 


2.0 


.62 


.967 


1.25 


0.34 


101.5 


65.0 


Table 


VII. 


Recorded 


and predicted limit cycle 


values for 


■ ramp 


inputs. 
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describing function predicts the ratio of the fundamental sinusoidal 
components of the input and output signals. In many of the runs con- 
ducted (large and small ratios of ) the input to the non-linear 

iA v 

element uas decidely non-sinusc idrl. This can be appreciated by 
examining the Brush recordings. It should be noted that the value 
recorded for the experimental value of X was simply ■§• the peak to peak 
amplitude of the relay input, while the value of X predicted by the 
describing function is the amplitude of the fundamental frequency 
comoonent. Aopendix II shows that for values of -JL- near 0,1 or 0.9, 

iXy 

and for values of fi near 90°, the pull-in angle approaches ISO 0 , thus 
requiring an erctrenely non-sinusoidal input signal with a peak to peak 
amplitude which approaches j the amplitude of the fundamental frequency 
component as illustrated by figure 4. 5. 

In section 4.2 and appendix II the limits for sinusoidal opera- 
tion are discussed and the area defining these limits is shown on the 
describing function curves in figure 4.2. It can be seen from the table 
of predicted and observed values that the errors are largest outside 
the area of sinusoidal operation, and that the predicted values become 
increasingly more accurate as the area of sinusoidal operation is ap- 
proached. 

Another consideration in evaluating the discrepancies between pre- 
dicted and observed values is the sensitivity of the describing function 
to hysteresis. The analog simulation of the relay employs two d.c. 
amplifiers in the analog computer, set to high gain values. It was not 
unco mon for the drift in these amplifiers to cause changes in the 
hysteresis setting of several tenths of a volt. In order to gain an 
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appreciation Tor :a ,uxtuue of errors introduced by s .all changes in 
hysteresis, a representative system containing fourth order linear com. 
ponents was examined,. The results of this examination are shoxm in 
table VIII and graphically aispla ed in figure 4.8, The extreme sensi- 
tivity of the amplitude to hysteresis changes is evident. It is also 
noteworthy that these errors are essentially invariant with the ratio 



B 




In conclusion, after careful analysis of the data recorded, it is 
apparent that the general validity of the proposed describing function is 
proven. Specifically, when intersection with a describing function curve 
occurs within the area defined for sinusoidal operation the amplitude 
and frequency of the limit cycle are predictable within 10m error, then 
the intersection occurs outside of this area proportionately greater 
errors are expected, but reasonable and useful, estimates can be obtained. 
4. 5 The Ba.no Describing Function for Two-sided Operation . 

If the sinusoidal component of the relay input is greater than 
the critical value of X defined by equation 4.12 or 4.19* the system 
will enter a limit cycle in which both sides of the relay are actuated 
during each cycle. In this case, the pull-in and drop-out angles must 
satisfy equation 4.20, which is obtained from the simultaneous solution 
of equations 2.20 and 1-12. 

4.22 +/S -o( 1 -/6 1 = ZyJL 

1 ' <v 2 

In order to evaluate the describing function for a relay with 

hysteresis it is necessary to obtain the simultaneous solution of equa- 
tion 4.22 and equations 1-3 through 1-11, This solution can be obtained 
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Table VIII, Errors in predicted amplitude and frequency resulting fro:.: 
snail relay hysteresis errors. Zero error condition is 
for p-d = 1,0, Data shown is for a system with fourth 
order linear components: 

2 

s(s +' V)' "(s + 2) (s" + ^i-) 
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ior constant values 



graphically by : lutiirg o( ~ 0 o( *,0 vs e 

of X, The value c vJ car. x ■ 1 be So-. insu fcr any combination of 

■q rr 

and X„ This Infer. ration can an turn be used to calculate de- 

YK 

v? 

scribing function cur ves which ari.ll be valid for taro-sided operation, 
Unfoi’tunately, this nethoa is ncc easily adapted to a non-dimension- 
alized form. Consequently it must normally be calculated individually 
for each relay to be investigated. 

The solution js slightly less complicated for a relay -without 
hysteresis. The transfer characteristics for such a relay are shown in 
figure 3. 1 „ In this case, °( - /S and o( = A,, Therefore the system 

of simultaneous equations which must be solved in order to evaluate W 
consists of only equation 4,23 and equations 1-9 and 1-11, 



4,23 0. - /3 X 




The equations for the coefficients of the fundamental components 
of the fourier series expansion of 0(t) reduce to 



4,24 a <5 — C 



4,25 b - flL ( C.:S 0 i m 

• rjj . £ 

The solution can be simplifies further by noting that 



4,26 cos $ cos 0 v = 2. cos ’?( 0 < + 0 ? )cosf-(^ - 0^) 



— 2 cosf( 0^+0 9 )cos 



2Yh, 



4,2? b, = iicos ;-( $ , x /3 0 )cos 

rr 



2YK. 



v 2 



?0 



Ail eve si -lmr oc • -ic o...ri je ^ rained for a system containing 
an ideal relay. In this case tile re is no dead zone, and ~ p 

Therefore equation 4.27 reduces to equation 4.23 and the describing 
function can be expressed, by equation 4.29. 
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It should be noted that there vail be a discontinuity between the 
describing function for one-sided operation and the describing function 
for two-sided operation whenever there is a significant hysteresis effect. 
Therefore, if the system is operating on only one side, and a small 
change in input or gain causes the error signal to change so that 
X > X„ , the system output may change radically. 

It is possible that this discontinuity may help to predict some 
types of "jump resonance" found in relay servos. Suppose, for instance, 
that the system is subjected tc a sinusoidal input, <9^ = As in ou^t, such 
that Au u t = .9 IXv2 c u , , I s aueh less than the natural frequencies 
of the system, the servo will behave as if the. input signal were a 
slowly changing ramp which has a maximum slope of Au> 0 at u» 0 t = 0 and 

decreases to zero am iot ~ ^ . 

° 2 

how suppose the system characteristics are such that the relay 
will operate on only one side when ^ 0 t = 0. Then the initial operating 
point will be determined by the ramp describing function for one-sided 
operation. The operating ooint mil tend to travel along the linear 
curve, with the instantaneous operating point being determined by setting 
B = Aui 0 cos u> 0 t. 
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If tho linear curve crosses the critical curve of for the 
relay used, the operating point will eventually reach a point where one- 
sided operation can no longer occur, and a two-sided limit cycle will 
result. During the next portion of the cycle, the servo will operate 
in accordance xjith the describing function for two-sided operation. The 
servo will return to one-sided operation when X decreases so that 
Du iL— (1 + sinc*p) , where both X and o< „ are determined by the 

Pn * 't 

tx.’c-sided describing function. 

In general, the magnitude of the system output will "jump" to 
a larger value when the relay shifts from one-sided operation to two- 
sided operation, and will "jump" to a smaller value whenever the relay 
returns to one-sided operation, due to the discontinuity between the two 
describing functions. There nay also be a change in the frequency of 
the limit cycle when the shift occurs. Both the frequency change and 
magnitude change can easily be predicted once the describing function 
curves have been plotted. 

4,6 Sxtension of the Rarro Describing Function to Other Types of I~on- 
linearities . 

The ramp describing function can be extended to apply to other 
types of non-linearities. In general, the fourier series expansion of 
the output wave will give several equations of the following types ; 

1 . A relationship for a 0 similar to equation 1-3. 

2. Relationships for a* and b^ similar to equations 1-5 and 1-7. 

3. Relationships for all angles involved, similar to equations 
1-3 through 1-11. 
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jjLnce e .aAuiv-ji c -jlies for any non-linearity, it can be 
combined with che foarier equivalent of a 0 to get a relation between 
the ramp magnitude and the angles involved, similar to equation 4.4, 

This in turn can be combined with the angle relationships to solve for 
the d„c. component, .4 Once b is evaluated, it is a relatively simple 
matter to evaluate the equations for a-j and b-j and substitute these 
values into equation 3=3 to obtain the describing function. 

It should be noted, however, that the transcendental equations 
involved will become more complex as the nature of the non-linearity 
becomes more complex, and the solution of these equations may become ex- 
tremely difficult. In addition, other types of non-linearities may have 
different modes of operation corresponding to the one-sided and two- 
sided modes of operation for the relay servomechanism, and there may be 
discontinuities between the describing functions for these different 
nodes. 

Thus a complete analysis will require the evaluation of a separate 
group of describing function curves for each node of operation, along 
with a set of limit curves to define the points at which the system wall 
shift from one mode of operation to another. 
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CHAP T2R V 

AHca:.c::DA;io::p for future i.:v33tigatio::s 



During the course of this investigation, several interesting 
problems were noted which appeared to be worthy of further study. These 
suggested subjects ares 

1 , A comprehensive study of the effect of different types of 
dissymmetry on systems containing different types of non- 
linearities , 

2. A more exhaustive analysis of the ramp describing function 
for a relay servomechanism in which the relay operates on 



both sides. This analysis should include an investigation 
of the discontinuity between the describing function for one- 
sided operation and the describing function for two-sided 
operation to determine whether or not these describing 



functions can be used to predict certain types of "jump 
resonance" effects, as suggested in section 4,5. 

3. A thorough study of the extension of the ramp describing 

function to non-linear devices other than the relay would also 
be valuable to determine the effect of ramp inputs on systems 
containing such elements. An outline for the procedure for 
such an investigation is given in section 4.6. 
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. APPENDIX I 



DERI VAT ICI'I OF THE GENERALIZED DESCRIBING FUNCTION 1 
FOR A DI3SYJ METRIC RELAY 



Input vs e output characteristics o 




Figure 1=1 
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_ . Fourier series expansion of output wave 
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Tor a nerfcct sine in rat, the angles can be defined as follows: 
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C. The describing fin': .cion. 



The generalized describing function is; 




If both sides of the relay are actuated on each cycle, equations 
1-4 through 1-7 can be used in the evaluation of equation 3«S without re- 
strictions. If the rale, is only partially actuated, the sane equations 
nay still be used grovi-'c • certain angles are defined as follows; 



If 
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n = 1,2 
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n = 1,2 
n = 1,2 
n = 1,2 



if 1-j = l r 2 = Y, equations 1-3, 1 - 5 * a nd 1-7 reduce to 

: ' 12 a o = W ( ' /3 2 + ^2 “ “ *1^ 

1-13 a 1 = ~(s in /# -| - sin ^ + sin /6* - sin ) 

1-14 = ^(cos j5 ■ j + cos c( ^ + cos /S 2 + cos o( 

If Y-j = Y£ = Y and the relay has no hysteresis, then = By 
and °( 2 = A 2 * ana equations 1-12, 1-13 and 1-1 4 reduce to 

1-15 a o = fr ( ^ 2 “ 

3*6 a — Q 

1-1 6 b-j - ~~(cos A 7 + cos ^ ) 

If Y.| = Y 2 = Y and the relay operates only on the positive side, 
then A 2 = ^ o ^ 7 < The subscripts can be dropped for convenience in 
writing, so that equations 1-12, 1-13* and 1-14 reduce to 



4 . 1 a 0 = | ( TT - /S - * ) 

4*2 a 1 = —(sin /S - sin o( ) 

1 rr 
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^ .LU. ±~. u.' nZ OF JI..730IDAL OPSRATIOii DURING 

Cr'iwt.’IC.: 0. Ox.L- 07 J FIDS 

.lien one relay operates on only one side, sinusoidal operation 



is possible 


only if o( < •■'0°. 
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For the case where > 2T , equation II-2 reduces to 
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Therefore if 3 — .yC:'^, sinusoidal operation is possible only if 
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